Effects of Disorder in FeSe : An A b Initio Study 



Prabhakar P. Singh 

E-mail: ppsingh@phy . iitb . ac . in 

Department of Physics, Indian Institute of Technology Bombay, Powai, Mumbai- 
400076, India 
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small amount of excess Fe substantially disorders the Fe-derived bands while 
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enhances the possibility of Fermi surface nesting, specially in FeSeo.sTeo.5, in spite 
of disordering the Se-derived bands, (iii) the electron doping through Co, Ni or Cu 
disorders the system and pushes down the Fe-derived bands, thereby destroying 
the possibility of Fermi surface nesting. A comparison of these results with the 
rigid-band, virtual-crystal and supercell approximations reveals the importance 
of describing disorder with the coherent-potential approximation. 
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1. Introduction 

The superconductivity in iron pnictides [U |2] and chalcogenides [3], with the 
interactions resulting from charge and spin degrees of freedom of electrons 
delicately balanced (31 [B], provide a unique opportunity to unravel the 
mystery of unconventional superconductivity [7J [51 [HI HDJ- From the experimental 
characterizations [11 of the superconducting state of the iron pnictides and 
chalcogenides as well as theoretical calculations [12_, it seems clear that the 
superconductivity in these alloys is not mediated by the phonons [T3]. Therefore, 
attempts are being made to understand the superconductivity in iron pnictides and 
chalcogenides in terms of spin-fluctuations and related theories |14j . 

In spin- fluctuation theories [5J[Tni[T3], in addition to the proximity of the system 
to a magnetic instability, its Fermi surface (FS) and, in particular, the FS nesting 
plays a crucial role in enhancing the magnetic interactions. In turn, by tuning the 
magnetic properties and the FS of the system by chemical substitutions, electron 
or hole doping, one can bring out the details of the superconducting properties of 
the system, and thereby close in on the exact nature of interaction responsible for 
superconductivity. Out of all the recently discovered, superconducting iron pnictides 
and chalcogenides, FeSe and its alloys with one of the simplest crystal structure are 
well-suited for such a study. 

In FeSe, the superconducting transition temperature T c of ~ 8 K [3 , increases up 
to 15 K with S and Te substitutions [TS] but decreases rapidly with electron doping 
using Co, Ni and Cu substitutions |151 [TBI 117] . The presence of excess Fe as well as 
the Se deficiency in FeSe is known to affect its superconducting properties |18) . In 
addition, with increase in pressure up to 8.9 GPa, the T c of FeSe increases to 36.7 K 

Previous theoretical attempts [22, 23, 24, 25, 26 at understanding the changes in 
the normal state electronic properties, as a prelude to understanding the changes in 
their superconducting properties, of several of the pnictides and chalcogenides upon 
alloying have used either virtual-crystal or the supercell or a series of ordered alloys 
or parametrized model Hamiltonian approach. The virtual-crystal and the supercell 
approaches are known to be quite inadequate to describe the effects of disorder in 
metallic alloys, especially in d-band metals |27| . In addition, an ab initio study of the 
effects of alloying in FeSe has been lacking so far. 

In an attempt to understand the changes in the normal state electronic properties 
of FeSe upon alloying, which may, in turn, provide some understanding of the changes 
in its superconducting properties, we have used Korringa-Kohn-Rostoker coherent- 
potential approximation method |27| in the atomic-sphere approximation (KKR-ASA 
CPA) [28 to carry out charge, self-consistent electronic structure calculations of FeSe 
and its alloys with S, Te, Co, Ni and Cu, including the presence of excess Fe and the 
deficiency of Se. We have also carried out a detailed comparison of these results with 
the rigid-band, virtual-crystal and supercell approximations. 

Based on our calculations, we find that (i) a small amount of excess Fe 
substantially disorders the Fe-derived bands near Fermi energy (Ep) while Se- 
deficiency affects mainly the Se-derived bands away from Ef, (ii) the substitution of 
S or Te on Se sub-lattice enhances the possibility of FS nesting in spite of disordering 
the Se-derived bands, specially in FeSe .5Te .5 alloy , (iii) the substitution of Co, Ni or 
Cu on the Fe sub-lattice disorders the system and pushes down the Fe-derived bands, 
thereby destroying the possibility of FS nesting. We also find that for describing 
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substitutional disorder the coherent-potential approximation is more reliable than the 
rigid-band, virtual-crystal or supercell approximation. 

2. Computational Method 

We have studied Fei +(5 Se with 5 = 0.0,0.01,0.03,0.06, FeSei_ x with x = 
0.01,0.03,0.06,0.12,0.18, FeSei_ y Te y with y = 0.25,0.5, FeSe . 9 S .i, Fe . 9 Co .iSe, 
Feo.gNig.iSe and Feo.giCuo.iSe in the tetragonal (P4/nmm) crystal structure, while 
Fei.oiSe has been studied in the orthorhombic structure as well. The substitution of 
10% Cu on Fe sub-lattice in Fei.oiSe results in the composition Feo.91Cuo.10Se, used 
in the cited experiment |17) and hence, in our calculations. To reduce the errors due 
to the ASA, we have introduced four empty spheres in the unit cell containing two 
Fe and two Se atoms. The atomic as well as the empty-sphere positions in the unit 
cell are: Fe (3/4,1/4,0; 2a), Se (1/4, 1/4, z Se ; 2c), El (1/4, 1/4, z = -z Se ; 2c) and E2 
(3/4,1/4,0.5; 2c), where El and E2 denote the two empty-sphere sub-lattices. The 
El empty sphere layer, in the same plane as the Se-layer, were used to accommodate 
the excess Fe in the calculations of Fei+^Se and Feo.91Cuo.1Se. To see the effects 
of incorporating the excess Fe in other interstitial regions, we have also used the E2 
empty-sphere site for the excess Fe in Fei 6Se and sites just below the Se-plane. 

To model the effects of disorder, we have used the CPA [29, 27J rather than a 
rigid-band, virtual-crystal or supercell method because CPA has been found to reliably 
describe the effects of disorder in metallic alloys [28, 27_. We used Barth and Hedin 
[30] exchange-correlation potential. The Brillouin zone (BZ) integration during self- 
consistency was carried out using a grid of 24x24x20 points in the BZ. The density of 
states (DOS) was calculated with a grid of 28x28x24 points in the BZ except for the 
pure FeSe and FeTe, where a grid of 36x36x32 points was used. For both DOS and 
spectral function calculations, we have added a small imaginary component of 1 mRy 
(1.5 mRy for FeSe and FeTe) to the energy. In the following figures the Se-derived s- 
band is not shown. In our calculations, the lattice parameters a and c were taken from 
experiments |151 [T51 [TBI IT7] while the theoretically relaxed z- values were taken from 
Ref. [T3]. For some of the intermediate alloys, we used the concentration- weighted 
average of the z- values. The lattice parameters of FeSe, FeTe and their alloys, used 
in the present calculations, are listed in Table [l] 

We have analyzed our results using the Bloch spectral function j^T] A(k,E), 
defined by 

A(k,E) = --%G(k,E), (1) 

where G(k, E) is the k-space coherent-potential Green's function, and k and E 
represent the wave vector and the energy, respectively, of the electron. The band 
structure along BZ symmetry directions was calculated by evaluating Eq. (JTJ for the 
given k points and the energy E in the given range. The Fermi surface in a given k- 
space plane was mapped by evaluating Eq. over a two-dimensional grid of 151 x 151 
k points in the plane at the Fermi energy E = Ep. The peaks in the spectral function 
A(k, Ep) form the Fermi surface of the alloy. All the plots of the band structures 
and the Fermi surfaces were carried out with the same value for the colormap in the 
range of to 30 (15 for the Fermi surface) with the minimum and maximum value 
represented by blue and red, and the intermediate values assigned colors varying from 
blue, light blue, green, yellow, orange and then red. Thus, a movement away from red 
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Table 1: The experimental lattice parameters a and c of FeSe, FeTe and their alloys 
used in the present calculations. 



Alloy 


a (A) 


c(A) 


Alloy 


a (A) 


c(A) 


FeSe 


3.765 


5.518 


Fe .9oCo .ioSe 


3.7637 


5.5043 


Fei.oiSe, FeSe .9 9 


3.7734 


5.5258 


Fe .9oNi .ioSe 


3.7713 


5.503 


Fei. 03 Se, FeSeo.97 


3.7787 


5.5208 


Fe .9iCu .ioSe 


3.807 


5.495 


Fei. 06 Se, FeSeo.94 


3.7747 


5.5229 


FeSeo.goSo.io 


3.763 


5.503 


Fei.i 25 Se 


3.7747 


5.5229 


FeSeo.75Teo.25 


3.7989 


5.9685 


FeSe .88, FeSe .875 


3.7693 


5.4861 


FeSe .5oTe .5o 


3.7909 


5.9570 


FeSe . 82 


3.7676 


5.4847 


FeTe 


3.8215 


6.2695 



and towards blue in the band structure and Fermi surface plots indicate reduction in 
the peak value of A(k, E). 

The formation energy, Ef orm , of excess Fe in Fei+^Se and Se vacancy in FeSei-^ 
alloys has been calculated as 

Eform — Ep ei+s Se ~ E FeSe - SE Fe (2) 
Eform = -^FeSei_ x — E FeSe + xE Se (3) 

where the superscripts e and v correspond to the formation energy of FeSe with excess 
Fe or with Se vacancy, respectively. The subscripts on the right hand side of Eqs. Q 
and |3]) denote the calculated total energies of the corresponding alloys except for E Fe 
and Ese which represent the corresponding free atom energy. 

The rigid-band calculations were carried out by first determining the amount of 
shift in Fermi energy required to either accommodate or remove a given number of 
electrons from the self-consistent potential of ordered FeSe. The Fermi surface in the 
rigid-band was then evaluated with the KKR-ASA CPA method at the shifted Fermi 
energy and using the self-consistent potential of ordered FeSe. 

For the virtual-crystal approximation |31| . we replaced the disordered sub-lattice 
by an ordered sub-lattice containing a virtual atom, and then carried out a charge, 
self-consistent calculation for the alloy with the virtual atom using the KKR-ASA 
CPA method. The virtual atom on a sub-lattice was constructed by adding the 
concentration-weighted average of the valence electrons of the respective atoms on 
the sub-lattice to the core electrons of the atom having higher concentration. For 
example, the atomic number of the virtual atom on the Fe sub-lattice in Feo.9Nio.1Se 
alloy is 26.2. 

In the supercell approximation for excess Fe in Fei.i25Se and Se vacancy in 
FeSeo.875, we used a 2a x 2a x c supercell of FeSe containing 32 sites (16 atomic 
and 16 empty-sphere sites) as shown in Fig. [I] The excess Fe in FegSes is placed at 
an empty-sphere site while vacancy in Fe§Se7 is formed by replacing a Se atom by an 
empty-sphere. The charge self-consistent calculations for the supercells were carried 
out using the KKR-ASA CPA method. The Brillouin zone (BZ) integration during 
self-consistency was carried out using a grid of 12x12x20 points in the BZ. The density 
of states (DOS) was calculated with a grid of 14x14x24 points in the BZ. For DOS 
and spectral function calculations, we have added a small imaginary component of 2 
mRy and 1 mRy, respectively, to the energy. We have also checked our results of the 
supercell calculations with the LMTO-ASA method. 



Effects of Disorder in FeSe : An Ab Initio Study 



5 




(a) (b) (c) 



Figure 1 : (Color online) The unit cell of (a) tetragonal FeSe and the 2x2x1 supercells 
used for (b) FegSeg and (c) Fe 8 Se7 calculations. The large spheres represent Fe (red) 
and Se atoms (yellow), and the smaller spheres represent empty-spheres (grey). The 
excess Fe in FegSes is placed at an empty-sphere site while vacancy in FesSe7 is formed 
by replacing a Se atom by an empty-sphere. 

3. Ordered FeSe and FeTe in the atomic- sphere approximation 

To be able to separate out the effects of disorder from the atomic-sphere- 
approximation-related changes, in this section we compare the electronic structure 
of ordered FeSe and FeTe alloys with the more accurate full-potential results. In Figs. 
[2]4] we show the band structure, the density of states and the Fermi surface of FeSe 
and FeTe alloys, which are similar to that of Ref. [13] except for the gap between 
Se- and Fe-derived bands around 2.5 eV below Ep in FeSe. Since it is already known 
that electronic structure calculations of FeSe are very sensitive to the height of the Se- 
layer above the Fe-plane, we have carried out full-potential, linear muffin-tin orbital 
(FP-LMTO) calculations for FeSe using the lattice parameters as above. We find 
that the gap predicted by FP-LMTO is «0.3 eV, which increases to «0.4 eV if the 
muffin-tin spheres are enlarged to atomic-spheres within the full-potential approach. 
In the atomic-sphere approximation, which uses additional approximations, including 
that of making the potential spherical, the gap is found to be «0.8 eV. However, the 
Fe-derived bands, which are responsible for superconductivity in these alloys, compare 
well with the more accurate FP-LMTO calculations as well as with that of Ref. [13] . 

Using the idea of FS nesting, it is possible to get some quantitative measure 
of the response of the system without evaluating the susceptibility. For example, if 
the FS around T point matches exactly with the FS around M point when displaced 
by a reciprocal space vector then the nesting is optimal. In the present context, the 
optimal nesting corresponds to FS at T and M points having matching radii and sharp 
Fermi surfaces (reflected by thin red lines in the figures). Any deviation, either from 
the matching radii or sharpness of the Fermi surfaces (reflected by diffused and/or 
broadened lines in the figures), generally, reduces the effect of nesting. 

The FS of FeSe, shown in Fig. [4] consists of two hole-like sheets around T point 
and two electron-like sheets around M point, both sheets being derived from xz (yz) 
and xy bands of Fe. In FeTe, there are three hole-like sheets around T point and only 
two sheets around M point as shown in Fig. [4] The Fermi surfaces shown in Fig. 
[2] reveal enhanced FS nesting at T-X-M plane in FeSe than in FeTe. Note that the 
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Figure 2: (Color online) (a) The band structure of tetragonal FeSe (top) and FeTe 
(bottom) along the BZ symmetry directions using experimental lattice parameters and 
theoretically relaxed zg e . The Fermi energy is indicated by the horizontal line at 
cV. (b) The colormap used for all the band structure and the Fermi surface plots as 
described in the text. 

two bands around M point in FeSe are not resolved in Fig. [2] In the following, we 
will see if the changes in the shape and the nesting of the FS of FeSe induced upon 
alloying can be used to understand the changes in the superconducting properties of 
these alloys within the framework of the spin-fluctuation theories. 

4. Effects of theoretical vs. experimental zs e 

The role of Se height above the Fe plane, given by zs e , plays a crucial role in 
determining the electronic properties, especially the magnetic properties [32 , of FeSe 
and its alloys. It is not surprising that z$ e plays such an important role in deciding 
the electronic properties of FeSe because a change in z$ e directly impacts the Fe 
d-orbitals which, in turn, affect the Fe-derived (i-bands around Ep. Most of the 
theoretical work on FeSe have used the experimental values of the lattice parameters 
a and c, and theoretically relaxed value for zg e . To see what role does zg e play in the 
present context, we have studied FeSe and its alloys using the experimental value of 
zg e = 0.266 and the theoretically relaxed value of zg e = 0.2343 |13| . 

In Fig. [5] we show the band structure of FeSe calculated with the experimental 
as well as the theoretically relaxed value of zs e . In both cases, the experimental values 
of the lattice parameters a and c were used, and the exchange-correlation potential 
was parametrized using the generalized-gradient approximation of Perdew et al. |33] . 
The use of experimental z§ e . which is larger than the theoretical value, allows most 
of the bands, including the Se-derived bands, to move up with respect to Ep. Such 
a movement results in three bands crossing Ep along T-X, M-r and Z-R. The use of 
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Figure 3: (Color online) The total (black, solid) and the sub-lattice-resolved DOS 
of tetragonal FeSe (top) and FeTe (bottom) obtained using experimental lattice 
parameters and theoretically relaxed zs e - The contributions from the Fe (red, dash) 
and the Se/Te (blue, dot-dash) sub-lattices as well as from the empty-sphere sub- 
lattices El (green, dot) and E2 (violet, double-dash dot) are shown. The vertical, 
dashed line denotes the Fermi energy. Note that the total DOS in the figure 
corresponds to per two atoms. 




r xz k 

Figure 4: (Color online) The Fermi surface of FeSe (top) and FeTe (bottom) in T-X-M 
and Z-R-A planes obtained using experimental lattice parameters and theoretically 
relaxed zs e - 
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Figure 5: (Color online) The band structure of tetragonal FeSe along the BZ symmetry 
directions using experimental lattice parameters with zg e equal to its experimental 
(top) or theoretically relaxed (bottom) value. In both the calculations, the exchange- 
correlation potential was parametrized within the generalized-gradient approximation 
as described in the text. The Fermi energy is indicated by the horizontal line at eV. 



4rr 




-6 -4 -2 



E(eV) 

Figure 6: (Color online) The total DOS of tetragonal FeSe with the experimental 
lattice parameters, including experimental zg e (blue, dash), and with theoretical z$ e 
(red, dot) calculated within the generalized-gradient approximation. The total DOS 
of FeSe (black, solid) as obtained in Fig. [3] is also shown. 
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theoretical z$ e , which reduces the Se-height above the Fe plane, affects the individual 
Fe-derived d-bands around Ep differently due to the orientation of the <i-orbitals. In 
particular, we find that the 3z 2 — 1 derived band along T-X and T-Z is suppressed 
more than the other two bands, resulting in only two bands crossing Ep along T-X, 
M-r and Z-R for theoretically relaxed z$ e - 

The effects of using experimental or theoretical z$ e get amplified in the total 
DOS as shown in Fig. [6j The higher value of zg e shifts the DOS towards Ep as 
well as redistributes some of the electronic states between and -2 eV. The Fe-peak 
in the DOS is closer to Ep for experimental zg e than for the theoretical zs e - From 
Fig. [6] we also find that the DOS of FeSe obtained with the Barth-Hedin exchange- 
correlation potential is essentially identical to the DOS obtained with Perdew-Wang 
|34| exchange-correlation potential within the local-density approximation. 

We have also calculated the Fermi surface of FeSe with experimental and 
theoretical zs e - As discussed above, in the context of the band structure of FeSe, 
we find three bands around T and Z points for experimental zs e - 



5. Electronic structure of disordered Fei +x Se and FeSe\_ x alloys 

Generally, the synthesis of FeSe results in an alloy with either excess of Fe or 
deficiency in Se [18, 35) [T7]. Measurements on such alloys have shown that Fei+^Se 
is superconducting only when 6 < 0.01 and the underlying lattice is orthorhombic 
|16| . On the other hand, Se-deficient alloys FeSei_ 2; remain superconducting for 
x < 0.18 in the tetragonal structure [3j[T8]. In order to understand the contrasting 
superconducting properties of Fei + ^Se and FeSei_ x alloys, we have studied their 
normal state electronic structure as a function of S and x with < <5 < 0.06 and 
< x < 0.18. 

The changes in the band structure, the density of states and the FS in Fei+aSe 
and FeSei_ x alloys with increasing S and x are shown in Figs. 7fl0 In the unit cell, 



the excess Fe atom, denoted by Fe(2) (The Fe at the stoichiometric site is denoted by 
Fe(l)), is kept in the Se plane with zp e — — zse, which is one of the two in-equivalent 
sub-lattices containing empty spheres. Our main conclusions are independent of the 
possible locations of the excess Fe atoms within the FeSe lattice, as was found in Ref. 
|18| . The increase in disorder induced by the presence of excess Fe can be clearly 
seen, as indicated by the diffused intensity, in the Fe-derived bands in Fig. |7ja). The 
Se-derived bands are left relatively untouched by the excess Fe. From the DOS of 
Fei+^Se alloys, shown in Figs. |8ja) and (c), and Fig. [9] we find that the excess Fe 
kept on the El sub-lattice creates states around Ep and around —2 cV. The possibility 
of the excess Fe becoming magnetic due to the increased local density of states at Ep 
can be clearly seen in the d-resolved DOS of Fe(l) and Fe(2) in Fig. [9] On the other 



hand, the increase in excess Fe quickly destroys the FS nesting as seen from Fig. 10 'a). 
Note that even for 5 = 0.01, the nesting is not optimal. 

For Fei.oiSe in the orthorhombic structure, the overall change in electronic 
structure with respect to the tetragonal structure is symmetry-induced and minimal. 
There are three bands crossing Ep around T point instead of two bands as is the case 
in the tetragonal structure, leading to some changes in the FS around T point. The 
DOS remains essentially unchanged with respect to the tetragonal case. 

The FeSei_ x alloys remain superconducting over a wide range of x~ values, 
indicating an unusually small impact on the bands around Ep due to vacancies in 



the alloy. Our results for FeSei_ x alloys, shown in Figs. 7fl0 seem to confirm it. The 
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(a) (b) 

Figure 7: (Color online) The band structure of (a) Fei+aSe and (b) FeSei^ with 6 and 
x= 0.01 (top), 0.03 (middle) and 0.06 (bottom) along the BZ symmetry directions. 
The Fermi energy is indicated by the horizontal line at eV. 



vacancies on the Se sub-lattice affect the Se-derived bands which are away from Ep, 
leaving the bands around Ep essentially untouched for small x. In addition to the 
lowering of Ep due to electron loss, for x > 0.06 the bands around Ep begin to get 
affected by disorder. For example, for x > 0.12, we find that the Se-derived bands 
have been substantially diffused. The loss of electron due to increasing Se vacancy 
moves Ep inward, seen clearly in Fig. |8jb), bringing the peak in DOS due to the 
stoichiometric Fe closer to Ep as shown in Fig. [8jd) , and thereby increasing the total 
DOS at Ep. The largest contribution to the DOS at Ep comes from the states of 
Fe(l) having x 2 — y 2 and xy symmetries. Such an increase in the DOS at Ep may 
lead to magnetism in the alloy [53] ■ We emphasize that Fe(2) atoms are responsible 
for possible magnetism in Fc 1+( 5Se, while Fe(l) atoms may lead to magnetic behavior 
in FeSei_ x . With increasing x, the gradual loss of possible FS nesting in Fei+^Se and 



FeSei-a; can also be seen in Fig. 10 If FS nesting is crucial for superconductivity in 
FeSei-z alloys then clearly the possibility of superconductivity is diminished for large 
x as the nesting is essentially destroyed as shown in Fig. |10| However, experiments 
|35II17| have shown that FeSe 1 _ :c for large x, contain impurity phases such as elemental 
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Figure 8: (Color online) The total and the sub-lattice-resolved DOS of (a) Fei+^Se 
and (b) FeSe!^ with S and x— 0.01 (black, solid), 0.03 (blue, dot) and 0.06 (red, 
dash). For FeSei-z, the DOS corresponding to x= 0.12 (magenta, double-dot dash) 
and 0.18 (green, double-dash dot) are also shown. The vertical, dashed line denotes 
the Fermi energy. The atom-resolved DOS at (c) the El-site of Fei+^Se and (d) the 
Se-site of FeSe^^ are also shown. 



Fe, leading to incorrect determination of Se content. Therefore, the Se content in 
superconducting FeSei_ a; alloys must be determined carefully. 

In both Fei+^Se and FeSei_ x , our results clearly show that the effects of excess 
Fe or Se deficiency on electronic structure of FeSe are very different, and that for 
x > 0.06 the FS nesting is essentially destroyed. 

6. Energetics of excess Fe and Se deficiency in FeSe 

To get an estimate of the energetics involved in the formation of FeSe alloys with 
either excess Fe or Se deficiency, we have calculated the formation energy of Fei+aSe 
and FeSei^a; alloys as described in Sec. [2j In Fei+^Se, we find that the formation 
energy Ej orm = -98, -94 and -87 meV/atom for S = 0.01, 0.03 and 0.06, respectively. 
The vacancy formation energy in FeSe 1 _ :E alloys is found to be E v form =115, 126, 
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Figure 9: (Color online) The d-resolved DOS at the Fe(l) (left panel) and Fe(2) (right 
panel) sites in Fei+^Se with 6 = 0.01 (black, solid), 0.03 (blue, dot), and 0.06 (red, 
dash). The vertical, dashed line denotes the Fermi energy. 




r x/. r r x/. k. 

(a) (b) 



Figure 10: (Color online) The Fermi surface of (a) Fei+aSe and (b) FeSei-a, with S 
and x= 0.01 (top), 0.03 (middle) and 0.06 (bottom) in T-X-M and Z-R-A planes. 
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Figure 11: (Color online) The total and the Fe(l) and Fe(2) DOS in Fei. 06 Se with 
theoretically relaxed zg e (left panel) and experimental z se (right panel). The DOS 
correspond to different positions of the excess Fe (Fe(2)) given by zp e — —z$ e (black, 
solid) with Fe(2) in the Se-plane, zp e = 0.7755 (blue, dot) with Fe(2) slightly below 
the Se-plane and zp e = 0.5 (red, dash) with Fe(2) in the middle of the two Se-planes. 
The vertical, dashed line denotes the Fermi energy. 



132, 134 and 136 meV/atom for x = 0.01, 0.03, 0.06, 0.12, and 0.18, respectively. As 
our calculations are done within the atomic-sphere approximation and without taking 
relaxation into account, these results are expected to change with the incorporation 
of full-potential and relaxation effects |23| . 

In order to check the effects of keeping the excess Fe in Fei+^Se at the various 
interstitial regions, we have studied Fei. 6Se alloy keeping the excess Fe (Fe(2)) in 
the Se-plane with zp e — —zsei slightly below the Se-plane with ZFe — 0.7755 and in 
the middle of the two Se-planes with zp e = 0.5 for theoretically-relaxed z$ e as well as 
experimental z se . As noted in Sec. [4] the DOS of FeSe and its alloys are expected to be 
different for theoretically-relaxed zs e and experimental z se , as is the case for Fei.oeSe 
in Fig. [TT] Not surprisingly, the DOS corresponding to zp e = — zg e and zp e = 0.7755 
are essentially identical for theoretically relaxed zg e . However, when the excess Fe 
is placed in between the two Se-planes instead of in the Se-plane, the resulting Fe(2) 
DOS's are somewhat different from each other, as shown in the bottom panel of Fig. 
11 These differences arise due to the relatively free space available to Fe(2) when 



zp e = 0.5. We also find that it is only with the inclusion of the Madelung potential in 
the Muffin-Tin-corrected total energy that Fe(2) prefers to go to zp e — —zg e instead 
of zp e = 0.5 by 11 meV/atom for experimental zg e . We must point out that a more 
reliable way of calculating the site-preference involving different structures is through 
the full-potential approach with relaxation included. 
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7. Electronic structure of disordered FeSeo.9So.i1 FeSe\- y Te y , Feo.9Coo.1Se, 
Feo.9Nio.1_Se and Feo.oiCuo.iSe alloys 

To understand the changes in the superconducting properties of FeSe upon 
substitution of impurities through the changes in the normal state electronic 
properties, we have studied the effects of substituting S and Te on Se sub-lattice 
and Co, Ni and Cu on Fe sub-lattice. In Figs. T2][T4 we show our results for the 



band structure, density of states and FS of FeSe .gSo.i, FeSeo.75Teo.25, FeSe .5Te .5, 
Feo.9Coo.1Se, Feo.gNi iSe and Feo.g1Cuo.1Se alloys. Based on the atomic size- 
mismatch, it is expected that the substitution of Te would lead to more disorder 
than the substitution of S on Se sub-lattice in FeSe. Indeed, we find that the effects 
of disorder are minimal due to S but substitution of Te creates states in the gap 
region with substantial disorder in the Se-derived bands. The states created in the 
gap between Se- and Fe-derived bands can be clearly seen in Fig. [l2{a) and (c). 
Surprisingly, for FeSeo.5Te .5 alloy, we find that the substitution of Te has rearranged 
the bands around Ep such that the possibility of FS nesting is maximized as shown 



in the bottom panel of Fig. 14 a). 

The presence of Co, Ni and Cu in Fe sub-lattice of FeSe is expected to affect the 
bands around Ep as well as move Ep up due to disorder and electron doping, thereby 
significantly changing its superconducting properties. As shown in Fig. |12[ b), a 10% 
addition of Ni or Cu moves Ep above the two bands around T point. The substitution 
of Cu seems to disorder both the Fe- and the Se-derived bands, as can be seen from 



the diffused intensity of the bands in Fig. 12 b). The addition of Co, Ni or Cu creates 



states in the gap region around -2 eV as seen from Fig. 13 b) and (d). The combined 



effect of disorder and electron addition diminishes the possibility of FS nesting in 
Feo.gCoo.iSe, and destroys it in Feo.gNio.iSe and Feo.g1Cuo.1Se as seen in Fig. [l4|b). 



8. Describing disorder : Rigid-band, Virtual-crystal, Supercell and 
Coherent-potential approximations 

In order to emphasize the importance of describing the effects of substitutional disorder 
accurately and reliably in metals, in general, and in FeSe and its alloys, in particular, 
we have studied selected FeSe alloys using the rigid-band approximation, the virtual- 
crystal approximation and the supercell approximation. In this section, we describe 
and compare our results, in terms of DOS and Fermi surface, of the rigid-band 
approximation, the virtual-crystal approximation and the supercell approximation 
with that of the coherent-potential approximation for the selected FeSe alloys. 

We have studied Fei+^Se and FeSei-^ with 5 and x= 0.01, 0.06 and 0.125, 
and Feo.gNio.iSe alloys in both the rigid-band approximation and the virtual-crystal 
approximation. For 5 and x — 0.125, studied for comparison with the supercell 
approximation, we have used the lattice parameters corresponding to S = 0.06 and 
x = 0.12, respectively We have also studied Fe g Se 8 and Fe 8 Se 7 as well as FegSeg 
in the supercell approximation. The selected alloys cover a wide range of systems 
with effects due to disorder involving addition of (/-electrons in Fei+^Se, removal of s- 
and p-electrons in FeSei-^, and similar atoms in Feo.9Nio.1Se. Generally, one expects 
the rigid-band approximation to work well for free-electron-like systems or systems 
consisting of similar atoms. On the other hand, virtual-crystal approximation is 
expected to work well for systems consisting of atoms from nearby columns of the 
periodic table such as Feo.gNio.iSe. The supercell approximation or its more efficient 
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Figure 12: (Color online) The band structure of (a) FeSeo.gSo.i (top), FeSeo.75Teo.25 
(middle), FeSeo.sTeo. 5 (bottom), and (b) Feo.9Coo.1Se (top), Feo.9Nio.1Se (middle) and 
Feo.91Cuo.1Se (bottom) along the BZ symmetry directions. The Fermi energy is 
indicated by the horizontal line at eV. 

version known as the special quasi-random structure |36| for describing the disordered 
alloy can be used to extract many properties of alloys such as the density of states 
and phase stability. 

In Fig. [15] we compare the total DOS of Fei +( 5Se and FeSei-^ with S and x= 0.01, 
0.06 and 0.125 calculated in the rigid-band and virtual-crystal approximations with 
that of the coherent-potential approximation. For S and x= 0.01, we find that the 
total DOS obtained in the rigid-band approximation is close to that of the coherent- 
potential approximation. For higher concentrations, the rigid-band DOS begins to 
differ substantially from the CPA DOS for Fei +< 5Se. Not surprisingly, in FeSei_ K the 
Se vacancy, involving the absence of s and p electrons, leads to a relatively small 
disagreement up to x — 0.06. 

For both Fei+^Se and FeSei-^, the DOS calculated with the virtual-crystal 
approximation differs substantially with the CPA DOS for 8 and x> 0.06, as can 
be seen from Fig. [l5j Note that for x— 0.125, the virtual-crystal approximation 
replaces the Se atom on the Se sub-lattice in FeSei-^ by a virtual atom with atomic 
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Figure 13: (Color online) The total and the sub-lattice-resolved DOS of (a) S and Tc 
substitutions on Se sub-lattice and (b) Co, Ni and Cu substitutions on Fe sub-lattice. 
The DOS in (a) correspond to FeSe .gSo.i (black, solid), FeSeo.75Teo.25 (blue, dot), 
FeSe .5Te .5 (red, dash), and in (b) Feo.9Coo.1Se (black, solid), Feo.9Nio.1Se (blue, dot) 
and Feo.91Cuo.1Se (red, dash). The vertical, dashed line denotes the Fermi energy. The 
atom-resolved, DOS of S and Te substitutions at the Se sub-lattice and Co, Ni and 
Cu substitutions at the Fe sub-lattice are shown in (c) and (d), respectively. 



number equal to 33.25. Thus, we find that the virtual-crystal approximation essentially 
fails to describe accurately the effects of disorder in the DOS of either Fei +( 5Se or 
FeSei_ 2 . for 5 and x > 0.06. However, the virtual-crystal approximation is expected 
to work better if the virtual atom is made out of atoms with similar atomic numbers 
as in Feo.9Nio.1Se. In Fig. [16] we compare the total and the Fe sub-lattice DOS in 
Feo.9Nio.1Se in the virtual-crystal approximation with the corresponding CPA DOS. 
Indeed, we find a good agreement in the DOS of Feo.9Nio.1Se. Some of the differences 
in the DOS of Fe .gNi .iSe around -1.9 eV and -1.0 eV highlight one of the shortcomings 
of the virtual-crystal approximation (as well as that of the supercell approximation) 
in assuming a periodic lattice with no energy-dependent electron scattering. The 
virtual-crystal approximation, as implemented here, cannot be applied to FeSeo.5Teo.5- 
A more appropriate approach for implementing the virtual-crystal approximation in 
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Figure 14: (Color online) The Fermi surface of (a) FeSeo.gSo.i (top), FeSeo.75Teo.25 
(middle), FeSeo.5Teo.5 (bottom), and (b) Feo.9Coo.1Se (top), Feo.9Nio.1Se (middle) and 
Feo.91Cuo.1Se (bottom) in T-X-M and Z-R-A planes. 

alloys may be through the use of psuedopotentials [37] . 

In Fig. [l5j the total DOS obtained for Fe 9 Se 8 and Fe 8 Se7, corresponding to the 
supercell approximation for 0.125 of excess Fe and 0.125 of Se vacancy per atom in 
FeSe, respectively, are also shown. The total DOS for FegSeg is in overall agreement 
with the corresponding CPA DOS of Fe! .125SC For Se vacancy, the total DOS of 
FegSe7 is very close to the fully-relaxed, full-potential supercell result of Ref. [55], as 
shown in their Fig. 2. A comparison of the supercell DOS of Fe§Se7 with the CPA 
DOS of FeSeo.875 shows some differences around Se-region and around Ep. Specially, 
the sharp peak close to Ep is missing in the CPA DOS. If the presence of the peak is 
confirmed by experiment, then it would indicate the need to go beyond the single-site 
CPA. The supercell approximation does provide the possibility of including several 
atomic environments, however, its inherent flaw of not including energy-dependent 
electron scattering limits its applicability. 

The k-resolved properties of alloys provide a more stringent test for the accuracy 
and the reliability of the various approximations for describing disorder that we have 
discussed so far. Therefore, we have calculated the Fermi surface of the selected 
alloys in the T-X-M and Z-R-A planes in the rigid-band, virtual-crystal and supercell 
approximations, and the results are shown in Figs. [T7][T9) The Fermi surface of 
Fei+^Se and FeSei_j; obtained in the rigid-band approximation, Fig. |17| agrees with 
the corresponding CPA Fermi surface, Fig. [lO] for 6 and x= 0.01 only. The Fermi 
surface obtained in the virtual-crystal approximation, shown in Fig. |18[ disagrees 
with both the rigid-band, Fig. [17] and the CPA, Fig. [10] results for 6 and x > 0.06. 
Similarly, for Feo.9Nio.1Se the Fermi surface in the rigid-band approximation is closer 
to the corresponding CPA Fermi surface than the Fermi surface obtained in the virtual- 
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Figure 15: (Color online) The total DOS of Fei+^Se (left panel) and FeSei_ x (right 
panel) with S and x— 0.01 (top), 0.06 (middle) and 0.125 (bottom) calculated in the 
rigid-band approximation (blue, dot), the virtual-crystal approximation (red, dash) 
and the coherent-potential approximation (black, solid). In the bottom panel, the 
total DOS obtained in the supercell approximation (green, double-dash dot) using 
FegSeg (left) and FegSey (right) are also shown. The vertical, dashed line denotes the 
Fermi energy. 
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Figure 16: (Color online) The total and the virtual-atom-site DOS of Fco.9Nio.1Se (red, 
dash), calculated in the virtual-crystal approximation as described in the text. The 
corresponding DOS in the rigid-band approximation (blue, dot), and the coherent- 
potential approximation (black, solid) are also shown. The vertical, dashed line 
denotes the Fermi energy. 
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Figure 17: (Color online) The Fermi surface, calculated in the rigid-band 
approximation as described in the text, of (a) Fei +( 5Se and (b) FeSei-^ with S and 
x= 0.01 (top), 0.06 (middle) and 0.125 (bottom) in T-X-M and Z-R-A planes. 
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Figure 18: (Color online) The Fermi surface, calculated in the virtual-crystal 
approximation as described in the text, of (a) Fei +( 5Se and (b) FeSei-a; with S and 
x= 0.01 (top), 0.06 (middle) and 0.125 (bottom) in T-X-M and Z-R-A planes. 
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Figure 19: (Color online) The Fermi surface of (a) Feo.9Nio.1Se, calculated in the 
rigid-band approximation (top) and the virtual-crystal approximation (bottom) in T- 
X-M and Z-R-A planes. In (b) the Fermi surfaces of FegSeg (top), FegSeg (middle) 
and Fe 8 Se7 (bottom) in T-X-M and Z-R-A planes, calculated using the supercell 
approximation, are shown. 



19 'a). Finally, we show the Fermi surface of 



crystal approximation as shown in Fig 
FesSeg, FegSes and FesSe7 in Fig. 19 b), obtained in the supercell approximation. 
Due to the folding in k x and k y directions, the Fermi surface of FegSe 8 and Fe 8 Se7 
cannot be compared directly with the CPA Fermi surface obtained using the primitive 
cell. Therefore, we also show in Fig. 19 b) , the Fermi surface of FegSeg. We find 



that the supercell approximation of FegSeg and FegSe7 predict a very different Fermi 
surface of disordered FeSe alloys than the coherent-potential approximation. 



9. Conclusions 



In conclusion, we have studied the electronic structure of Fei+^Se, FeSei-^, 
FeSei_ y Te y , as a function of S, x and y as well as FeSeo.gSo.i, Feo.gCoo.iSe, 
Feo.gNio.iSe and Feo.91Cuo.1Se alloys, our results show that (i) a small amount of 
excess Fe substantially disorders the Fe-derived bands near Ep while Se-deficiency 
affects mainly the Se-derived bands away from Ep , (ii) the substitution of S and Te 
enhances the possibility of FS nesting in spite of disordering the Se-derived bands, 
specially in FeSeo.5Teo.5 alloy , (iii) the substitution of Co, Ni or Cu disorders and 
pushes down the Fe-derived bands, thereby destroying the possibility of FS nesting. 
We also find that the coherent-potential approximation is more reliable than the rigid- 
band, virtual-crystal or supercell approximations for describing substitutional disorder 
in FeSe alloys. Thus, within the framework of spin-fluctuation theories, our results 
provide a consistent basis for understanding the superconducting properties of FeSe 
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alloys. 
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